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SUMMARY 

A comparison was made between the calculated design performance 
of a gas- turbine stator blade and its performance in a sector of an 
annular cascade tunnel. Information was obtained with regard to the 
three-dimensional effects that occur; the influence these effects have 
upon various performance parameters are presented. 

The gas velocities on the blade surfaces were computed by the 
stream- filament method and compared with the experimental values. 

The calculated values satisfactorily agree with the experimental 
values . 

The design conditions of free-vortex flow were not obtained. 

The experimental tangential velocities were greater than the design 
values at the tip of the blade but lower than the design values at the 
root of the blade. 


INTRODUCTION 

In the design of gas-turbine stator blades, compromises are 
generally made between aerodynamic considerations and such practical 
considerations as service life and cost of manufacture. In order for 
a turbine to perform efficiently and to meet design conditions of 
power and weight flow of gas, however, the blades must be designed by 
a procedure that accurately predicts blade performance. 

The blade profiles are generally designed according to two- 
dimensional-flcw theory. The velocity distribution along the blade 
surface must be such that flow separation or excessively high local 
velocities will be avoided. This velocity distribution may be 
determined by an incompressible potential- flow theory, which is appli- 
cable for any solidity or camber of the blades (references 1 to 4). 
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These methods, however, do not consider compressibility and are 
mathematically complex. Stream- filament theories that consider 
compressibility may be used to design blades of high solidity. 

These methods involve graphical integration of the flow equations 
(reference 5). Satisfactory experimental verification of either 
the potential-flow theories or of stream-filament theories have 
previously been lacking. 

The design performance of turbine stator blades based on free- 
vortex flow as measured in a sector of an annular cascade tunnel 
was determined at the NACA Lewis laboratory. The blade- surface 
velocities were calculated by an approximation, which was developed 
at the NACA Lewis laboratory, of the stream-filament method of ref- 
erence 5 and are compared with experimental values. 


SYMBOLS 

The following symbols are used in this report: 

C curvature, (l/ft, unless otherwise noted) 

Cp specific heat at constant pressure, (ft-lb/(slug) (°B) ) 

c v specific heat at constant volume, (ft-lb/(slug) (°B) ) 

exp base of Napierian logarithmic system e raised to 

power in parentheses following exp 

n distance from blade suction surface to any point on 

velocity-potential line, (ft, unless otherwise noted) 

n,-, length of velocity potential line from suction surface 

to pressure surface, (ft, unless otherwise noted) 

P pressure, (lb/sq ft absolute) 

r radius, (ft, unless otherwise noted) 

S blade- surface length, (ft, unless otherwise noted) 

T temperature, (°E) 

V veloc ity , ( ft/ sec ) 

V cr critical velocity, 2c p T^, (ft/sec) 
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W 


a 

r 

7 

v 

P 

T 

* 


weight of gas flowing across velocity-potential line 
per unit depth, (slug/sec) 

v = fzzi _v_ 

V cr 

angle between normal to cascade axis and flow direction 
of gas, (deg) 

circulation, (sq ft/sec) 

ratio of specific heats, c p/ c v 

local deflection angle, (deg) 

density, (slugs/cu ft) 

blade pitch, (ft, unless otherwise noted) 
velocity potential, (sq ft/sec) 
stream function (slugs/sec) 


Subscripts: 

0 ambient static conditions 

1 suction (convex) surface 

2 pressure (concave) surface 

e cascade exit 

i cascade entrance 

m value at point of average streamline curvature 

r component in radial direction 

s static conditions 

t stagnation conditions 

u component in tangential direction 

x component in axial direction 
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Parameters : 
a 

b 

f 

F(t x ) 

F(t 2 ) 

g 

J 

K 

t 

y 

u 


n o c i 


£1 

AC 


- Z m> 


2-7 

7-1 


/2-y > 

\t -ly 


exp 


(t 2 )dt 


exp 


(t 2 )dt 


3 /2 2~7 

# (1-z.P 


exp 


n 


- m (c2 - O 


dC 

AC 


exp 


w* - tf) 


dC 

AC 


C . - 


2AC 


n 0 AC 


2 

W 


n /2c T, 
o a/ p t 
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Ps V x 


p t 

V cr 


dP 

r 

s 

p t 

dr 

p,. 

/V \ 2 

s 



w 

r 

p s av r 


v cr P t dt 


weight-flow parameter 


radial pressure-gradient parameter 


centrifugal-force parameter 


radial-acceleration parameter 


BLADE DESIGN 

The turbine stator considered was designed for the following 
conditions: inlet-gas pressure, 8200 pounds per square foot; inlet- 

gas temperature, 1900° R, gas flow, 40.4 pounds per second; power 
output, 4300 horsepower, and rotor-blade speed at the 10.35-inch 
radius, 1000 feet per second. The weight-flow parameter, critical 
velocity ratio, and gas-flow angle at three radial stations are 
listed in the following table: 


Station 

Radius 

(in.) 

Weight -flow 
parameter 

Critical velocity 
ratio 

fei 

Gas discharge 
angle 

a e 

(deg) 

(Ps V x \ 

/Ps V x \ 

Wor). 

Vt V crj e 

Root 

9.3 

0.212 

0.233 

0.878 

67.8 

Pitch 

10.35 

.212 

.249 

.800 

65.6 

Tip 

11.5 

.212 

.260 

.737 

63.3 


These values were computed on the assumptions of radial equi- 
librium, constant gas energy and entropy, and free-vortex flow from 
root to tip. These assumptions result in theoretically constant 
axial velocity at all radii (reference 6). Flow-area blockage due 
to boundary- layer development at the blade root and tip was con- 
sidered by assuming, in accordance with the experimental data of 
reference 6, a boundary-layer-displacement thickness of 0.^ and 
0.1 inch at the inner and outer shroud, respectively. 
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In establishing a velocity distribution over the blade sur- 
face, flow acceleration over as much of the blade surface as 
possible is considered desirable. Thus, blade shapes at three 
radial stations, root, pitch, and tip (table I) were designed with 
the stream-filament theory developed in appendix A by the procedure 
outlined in appendix B. Values of local deflection angle at cas- 
cade entrance, and local deflection angle at cascade exit 

of 5° and 2°, respectively, were used at all radii. The values 
were graphically obtained from previous designs based on potential- 
flow theory ( reference 3 ) . 

The blade coordinates and other dimensions are shown in table I. 
Blade shapes for a typical blade section with the velocity-potential 
lines used in calculating surface velocities are shown in figure 1. 
The calculated surface velocities are plotted against blade-surface 
length in figure 2. 


APPARATUS AND PROCEDURE 

The investigation was conducted with five blades mounted in a 
sector of an annular cascade to simulate three-dimensional flow. A 
sketch of the experimental equipment is shown in figure 3. 

The experimental conditions were different from the design 
conditions in that the blades were designed for an inlet-gas tem- 
perature of 1900° R and the experiments were conducted at an inlet- 
gas temperature of 580° R. The Reynolds number (based on the pitch 
of the blade) for the turbine blades at design conditions was 
600,000; whereas, for the cascade experiments, the Reynolds number 
was 800,000. The effect of this difference was considered small. 
The blades were, however, investigated at the design value of 
critical velocity ratio V/V cr near the pitch section of the blade. 

The pressure distribution about the center blade of the cas- 
cade was measured with 36 static-pressure tubes of 0.030- inch 
diameter located as indicated in table II. Pressures were measured 
with mercury manometers. 

A survey probe, which is shown in the enlarged portion of 
figure 3, was used to determine the angle at which the gas dis- 
charged from the cascade of blades. The gas angle, measured from 
the axial direction, was read on a protractor. The center tube of 
the probe was used to indicate the stagnation pressure; the static 
pressure was read by static-pressure holes located 9 diameters 
downstream of the tip of the probe. The position of the end of the 
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probe was varied in the discharge field by manipulation of the 
traversing mechanism. The static pressures at the inner and outer 
shrouds were read by wall taps located midway between blades in 
line with the trailing edge. 

A preliminary survey of the tunnel was conducted to assure a 
uniform pressure field upstream of the blades. The stagnation- 
and static-pressure distribution across the tunnel along the sur- 
vey paths (indicated in fig. 4) is shown in figure 5 for a plane 
1.5 chord lengths upstream of the blades. Inasmuch as no boundary- 
layer- control slots were provided, the boundary-layer development 
upstream of the blades along the inner and outer shrouds was inves- 
tigated. The approach was long enough to provide uniform flow and 
short enough to prevent excessive boundary- layer growth, as seen 
in figure 5. The variation of weight-flow parameter at the cascade 
entrance with the ratio of stagnation pressure at the cascade 
entrance to ambient static pressure is shown in figure 6. The test 
condition is indicated. 

The following measurements were made to determine blade 
perf ormance : 

1. Stagnation and static-pressure surveys, both upstream and 
downstream of blades 

2. Stagnation upstream temperature (held constant at 580° E) 

3. Static -pres sure distribution about center blade of cascade 

4. Survey of angle of discharge of blade 

The downstream surveys were made in a plane 0.1 chord length 
behind the blades . The surveys were made close to the blades 
because the inner and outer shrouds extended only 1 chord length 
downstream of the blades. 


I 


RESULTS AND DISCUSSION 

The stagnation and static pressures at simulated design con- 
ditions in a survey plane 0.1 chord length downstream of cascade at 
points midway between blades at a number of radial stations are 
shown in figure 7. The construction of the survey equipment limited 
the positioning of the probe to radii 0.3 inch less than that of 
the outer shroud. The data of figure 7 are plotted as the ratio 
of stagnation- to-static pressure in figure 8. The design value of 
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pressure ratio at any radius is the pressure ratio that would be 
obtained with radial equilibrium of pressures and free-vortex flow 
at design conditions. The experimental results show a pressure 
ratio at the root 2 percent lower than the design value and a 
pressure ratio at the tip 4 percent greater than the design value. 

Discharge angle. - The gas-discharge angle measured from the 
axial direction at a number of radial stations in a plane 0.1 chord 
downstream of the blades is shown in figure 9(a). In general, the 
air was turned approximately 1° more than the design value. The 
gas was considerably overturned in the boundary layer near the 
inner shroud. This region of severe overturn is followed by an 
adjacent region at a larger radius in which the gases were turned 
to the design leaving angle. Near the outer shroud the gases were 
slightly underturned. 

Although no data were taken at the outer shroud, the gases 
were probably overturned in this region. The overturning near the 
inner shroud is probably caused by secondary flow from the high- 
pressure region on the pressure surface of the blade to the low- 
pressure region on the suction surface as explained in reference 6. 

The variation in discharge angle in the circumferential direc- 
tion at the 10.35-inch radius is shown in figure 9(b). The gases 
near the suction surface of a blade were turned about 8° more than 
the gases near the pressure surface of the adjacent blade. Inas- 
much as the included angle between suction and pressure surfaces 
at the trailing edge is 10°, some difference between discharge 
angle near the two surfaces was to be expected. 

Vane-surface-velocity distribution. - The experimental- and 
calculated-design-velocity distribution about the root, pitch, and 
tip sections of the blade are shown in figure 2. The experimental 
values of velocity on the blade surface were computed from experi- 
mental values of static pressure, assuming no loss in total pres- 
sure in the region outside the boundary layer on the blade surface. 
With this assumption, the surface velocity at a point on the blade 
surface is related to the static pressure at the same point as 
follows : 



where the total pressure is measured upstream of the blade row. 
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In general, the experimental velocities are higher than design 
values on the pressure surface and lower on the suction surface. 

The differences sire more pronounced at the root and tip sections 
than at the pitch section; however, these sections are at the edge 
of the boundary layers on the inner and outer shrouds. The shear 
stresses in the boundary layer cause an increase in entropy and 
deviation from the theoretical velocities. 


Tangential- and axial-velocity distribution. - The gas -discharge 
velocity components were computed from the pressure ratios of fig- 
ure 8 and the discharge angles of figure 9(a), assuming that the 
radial component of velocity is small enough to be neglected. 

The experimental and design values of critical velocity 
ratio V/V cr are plotted against radius in figure 10(a). The 
measured values are 3 percent higher than design near the blade tip 
and 3 percent lower than design near the blade root. 


The measured values of tangential component of critical veloc- 
ity ratio V u /V cr are 3 percent higher than design near the blade 

tip and 3 percent lower near the blade root, as shown in figure 10 (h) 
Vortex distribution of V u /V cr was not obtained; but inasmuch as 

the values are, in general, greater than design values, the turbine 
power will probably exceed the design power by about 2 percent if 
the rotor is assumed to perform as designed. 


Constant axial velocity was not achieved (fig. 11(a)); the 
experimental values were about 6 percent lower than the design 
value at the 10- inch radius and 6 percent higher near the blade 
tip. 


The weight-flow parameter 


P V 
M s x 


obtained from the data is 


M t cr 

compared with the design values in figure 11(b). This parameter 
is quite close to the design value at the blade sections at which 
the blade- surface-pressure distributions were obtained; at the 
10-inch radius, however, this parameter is 6 percent lower than 
the design value. Although no accurate computation of the experi- 
mental boundary- layer thicknesses at the blade root and tip can be 
made because of insufficient data describing conditions near the 
tunnel shrouds, the check between design weight flow and experi- 
mental weight flow is considered close enough to say that the 
design allowance for boundary- layer growth at the shroud was 
adequate. The principal reason that vortex flow was not achieved 
is because radial equilibrium of pressure was not attained. The 
heavy boundary layer at the inner shroud also tends to deflect the 
air toward the outer shroud and distorts the flow near the inner 
shroud. 
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Radial-equilibrium consideration. - As stated in the discus- 
sion of blade design, the blades we re designed on the assumption 
of radial equilibrium of pressure. This condition may be stated 
as follows: 


d£s . P S V u 2 
dr r 

This equation states that the centrifugal force on the air particles 
is balanced by a pressure difference across the particles. However, 
because radial flow must take place as the gases pass through the 
passage between blades if the exit axial velocity is to remain 
constant at all radii, the gas will be accelerated radially outward. 
These radial acceleration forces can be considered by writing 

£« . Ps T U 2 . . 
dr r s dt 

This expression may be written in dimensionless form for convenience 
as 

7+1 r ^s _ £s \ r £s ^r 
2 7 P t dr p t \v cr y ■ v 2 p t dt 

cr 


The design assumption in effect neglected the radial acceleration 

dV_ 


of velocity term involving 


dt * 


The magnitude of the radial-acceleration parameter involving 
the rate of change of radial velocity with time at the experimental 
conditions is shown in figure 12. The radial -pressure-gradient 
r 0P S 

parmeter __ vbs calculated by graphically differentiating ^ 

the lower curve of figure 7. The centrifugal-force parameter — 1 = — , 

p t Ycr / 

was computed from figures 10(b) and 7 by assuming isentropic flow. 


The radial-acceleration parameter 


r Ps 

V 2 p. dt 
cr z 

Pc /V \2 


was then eval- 


uated by subtracting Z±L JL from -S. _JL 


27 dr 


Pt 


Although 


' cr / 
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av dp 

Is nearly as large as - — ~ 


2 dt 


'cr 


%\ 2 


was small (about 


effect on — 
p t 

tion from free- vortex flow was small. 


cr/ 


near the blade tip, the 
4 percent) and the devia- 


SUMMARY OF RESULTS 

From a comparison of the calculated design performance and 
the experimental performance of a gas-turbine stator blade, the 
following results were obtained: 

1. The stream-filament theory is adequate to determine gas 
velocities on the surface of high-solidity blades. For the blades 
investigated, the calculated surface velocities satisfactorily 
agreed with the experimental values. 

2. The experimental value of tangential component of critical 
velocity ratio was 3 percent greater than the design value at the 
blade tip and 3 percent less than the design value at the blade 
root. These deviations were caused by lack of radial equilibrium 
of pressures at 0.1 chord downstream of the turbine blades and the 
thick boundary layer at the inner shroud. 

3. Constant axial velocity at all radii was not achieved. 

The experimental axial velocity was 6 percent higher than the 
design axial velocity near the blade tip. 

4. Although the radial-acceleration forces were large, the 
deviation from free- vortex flow was small. 


Lewis Flight Propulsion Laboratory, 
National Advisory Committee for 
Cleveland, Ohio, September 


Aeronautics, 
13, 1948. 
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APPENDIX A 


STREAM-FILAMENT THEORY 

The tvo- dimensional flow between curved surfaces may be 
studied by the stream- filament theory if the channel length is 
great compared to its width. This method consists in combining 
the conditions of irrotationality of flow with the principle of 
continuity of flow. 

Determination of velocity variation across channel. - In the 
general case of the motion of a compressible, frictj onless, perfect 
gas between curved surfaces (fig. 13), at any point in the channel 
the fluid rotates about an instantaneous center at radius r with 
velocity V. Inasmuch as the fluid is considered irrotational, 

= _ X = - CV (1) 

dn r 


The derivative is negative, inasmuch as the positive direction 
is chosen as that from the suction surface to the pressure surface. 

The curvature of the streamline is assumed to vary linearly 
from the suction surface to the pressure surface and is considered 
positive if the center of curvature is below the surface. 


C = c x 



( 2 ) 


Then 




dC 



(3) 


(4) 


If equation (1) is rewritten in terms of curvature and equation (4) 
is used. 



C dC 


d V 
V 


(5) 


1026 


1026 


NACA TN No. 1810 


13 


When equation (5) is integrated. 


log e V = - 


HoC2 


+ constant 


( 6 ) 


2AC 


or 



(7) 


where V and V m represent the velocity at any two points along 
the potential line. 

Evaluation of velocity from mass flow. - If the velocity at 
any point along a potential line is known, then the rest of the 
velocities along the line can be determined, as indicated in 
equation (7). From the use of the continuity equation, the veloc- 
ity at one point may be determined. 

The differential quantity of mass flow per unit depth of 
channel 


dW = p s V dn 


( 8 ) 


or in terms of density and temperature at stagnation 


1 



(9) 


Now, let 



( 10 ) 


then. 
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® - (1 - z) 5 " 1 Jz 

i i 

(1 - z ) 7 ' 1 = (1 - Z m - z + z m ) 7 " 1 


( 11 ) 


= (1 



Z - Zmf' 1 
" 1 ' V 


(12) 


where 





2 c T + 
P t 


(13) 


The second term in equation (12) may be expanded by the binomial 
theorem, choosing A /Z^ at the point of average streamline curva- 
ture (channel center) so that Z - will be small, and neglect- 
ing all but the first two terms of the expansion, then 


1 

7-1 

(1 - Z) = (1 



1 + _1_ 1 

f Zm N 

| . _i_ 1 

f z \ 

L 1 7-1 1 

VI - 

1 7-1 1 

\1 - V 


Let 


= (1 



2-7 

7-1 


1 + ^ 



Z 

7-1 


f = (1 



(14) 


(15) 


and 
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Then, 



Then equation (ll) becomes 


dW 

p t n oV 2 c p T t 



Let 


d = 


p t n o V 2 °P Tt 


dW 


_ W 

Pt n o J 2 c p T t 

If equation (18) is integrated. 



but 


(16) 


(17) 


(18) 


(19) 


( 20 ) 


fT _ X 


and from equations (7) and (4) 
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let 


and 




exp 





Then equation (20) becomes 


p = f J - gK 


( 21 ) 


( 22 ) 


(23) 


or 

f = S f + E ( 24 ) 

Equations (22) and (23) may be evaluated for J and K for C^Cg 
by letting 


Then, 


and 


n o C 2 
2AC 



dC 

AC 




(25) 
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Then, 


and 


J = 



exp I* dt 


exp (t^) dt 


(26) 



The terms t^ and are known from the channel dimensions 

PX 

and a table of / exp(t^)dt can be found in reference 7. 

Jo 

Evaluation of velocity at channel center. - The factors f 
and g are functions only of ^2^ and 7 (equations (15) and (16)). 

Thus, for given values of Jz^ and 7 , curves of f plotted 
against and g can be made, as shown in figure 14. In order 

to evaluate a line is drawn at tan”' 1 ' ^ intersecting the 

ordinate of the plot of f against g at (fig. 14). The line 

will, in general, intersect the curve at two points, one inter- 
section corresponding to the subsonic solution and the other, the 
supersonic solution. If the line is tangent to the curve, the 
channel is choked. If the line lies above the curve, no solution 
exists inasmuch as p is larger than physically possible. The 
value of is then found by drawing a horizontal line from the 

subsonic or supersonic intersection point, whichever is desired, to 
the curve of f against then is read on the abscissa 

as shown in figure 14. 
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The incompressible value of /Jz^ is determined by setting 
g = 0 that is, 

5 -'-^ ( 5 <0 - 2s ) 


Then 


Pt n o J ’ Vn 

Inasmuch as the flow is considered incompressible 


Pt = Ps 


V m = 


W 


p_ n J 
s o 


Evaluation of J and K. - A chart for the evaluation of J 
and K/j to be used in finding /Zm would be a simplification 
over evaluation from tables. A chart can be constructed where J 
and K/J are functions of the surface curvature and channel width, 
inasmuch as from equation (26) 


Let 


Then, 


J 1 = 


exp (- tjn 2 ) 


AC Dq 


F (t x ) - F (t 2 jj 



= 


AC 


C 2 + Ci 


(28) 


(29) 


(30) 


and 


F(t x ) = F (t m + |) 


(31) 


% 
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Now, let 


Then, 


“o C 1 

a = t-L y = — 2 ~ 


b y AC 


i (2b - l) 2 


and equation (34) becomes 


J = 1 + 


12b 


1 + h (2b - 1) ‘ 


F(t 2 ) = T - |) (32) 

These functions may be expanded in a Taylor's series. 

- 2 [© »■ (W - (I ) 3 t 51 ♦ • ■ ■ • ♦ (f f ^ 

(33) 

where n is an odd positive integer and 

F '( t m ) = ex P (tm 2 ) 

F ' ' ' ( tjjj) = ( 2 + 4 tjjj^ ) exp ( tjjj^ ) 

By use of the first five terms of the Taylor's series, 

J = 1 + (1 + 2 V 2 ) + ^ (3 + 12 t,* +4^) (34) 


(35) 

(36) 


1 ( a 


480 Vb 


3 + 3(f) (2b-l) 2 + |(f) 2 (2b - 1)' 


(37) 


and by the same method but using seven terms of the Taylor's series 
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K = 1 + 


4b 


I (!) <» - «' 


+ 1 > + 


160 


(iff! [(I) (**-»*] 


+ 3 (i) (2b - .^5 (|) I f(|) (2b - 1) 


f [l (2b ' 1)2 ] 2 * f(|> < 2b - 1 > 2 + 5 


(38) 


Equations (37) and (38) apply for negative or positive values 
of b, that is, for C^> Cg or Ci< Cg but not for = Cg . 

If Cl = C2, the expressions for J and K are 


J = 


slnh a 
a 


sinh 5a 
* 3a 


t 


These equations can easily be verified by integrating equa- 
tion 1 vith the curvature C constant and substituting in equa- 
tions (22) and (23). 

Charts of J and K/j are shown in figures 15 and 16, respec- 
tively, as a function of 


and 



Evaluation of channel -surf ace velocities, - The relation 
between the velocity at the point of the center of the channel and 
the velocities along the surface of the channel (figs. 17 and 18) 
can be derived, inasmuch as from equation (21) 


% 
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(39) 

(40) 


The velocities at the channel wall, that is, at the ends of 
the velocity-potential line considered, are computed from the value 
of Zjjj by use of equations (39) and (40). 


The simple stream -filament theory presented may be used to 
compute the velocity distribution on the surface of any channel to 
which the assumptions used apply. 
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APPENDIX B 


BLADE-DESIGN PROCEDURE 

The velocity distribution for cascades of high solidity may 
be determined by use of the stream-filament theory presented in 
appendix A. Blades may be drawn by inspection to fit the velocity 
vectors specified by the vector diagrams. Starting the blade 
design by first drawing the camber line seems to be practically 
desirable. The ends of the camber line are drawn tangent to a 
line at an angle of a e + v e from the horizontal at the trailing 
edge and to a line at an angle at the blade entrance 

(fig. 19) . The camber line may be symmetrically drawn about 
line X, vertically drawn through the camber line where its slope 
is zero. 

A symmetrical airfoil is superimposed upon the camber line, 
thus a first approximation of the blade shape is formed. The maxi- 
mum thickness of the symmetrical airfoil is placed at line X 
(fig. 20). A circle forming the leading edge is drawn with a 
diameter equal to about 30 percent of the maximum thickness of the 
foil. The trailing edge is drawn with an included angle of about 
12°. The maximum thickness of the foil will, of course, depend on 
the toted, tunning angle + a Q and the desired surface velocities. 

An orthogonal network of streamlines and velocity-potential 
lines are drawn in the channel formed by two airfoils. The network 
may be extended ahead of and behind the cascade so that the surface 
velocities may be computed at any point on the blade (fig. 19). 

The computed surface velocities in regions from to the rear 

stagnation point and from 0g to the front stagnation point are 
approximate, inasmuch as the channel is bounded by an assumed 
streamline in these regions. The velocities at both ends of each 
velocity-potential line are computed by use of the methods of 
appendix A. 

ja ru 

In general, the value of — - — at the channel entrance will 

not be equal to that at the exit because of deviations from two- 

dimensional flow. For convenience, this term can be considered to 

vary linearly with the axial depth of the blade, as shown in 

figure 21. The following relations, based on the general energy 

equation and the isentropic relations for a perfect gas are used 

u n 
H o 


to evaluate 


T 
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1 


7-1 

7+1 






7-1 

7+1 




(41) 


(42) 


With the flow network of velocity potentials and streamlines 
drawn and the value of p known for each velocity potential line, 

the surface velocity may be computed. The values of f-X_\ and 

V cr /j, 

are plotted against blade- surface length and S 2 , 

;cr/2 

respectively. (See fig. 22.) The velocity potential f is then 
determined by graphically integrating the velocity-distribution 
profile. 


0 - 


V ds 


(43) 


The velocity potential is plotted against blade-surface 
length S (fig. 22). Inasmuch as the absolute value of 0 is 
not important, the value of 0^, the velocity-potential line con- 
necting and Sg at the rear stagnation point on Sg, is set 

equal to zero on both and Sg. This plot provides a convenient 

method for relocating the velocity-potential lines, inasmuch as the 
value of S-j^ and S 2 for constant values of 0 are easily found. 

The circulation T around the airfoil is specified by the 
vector diagram and is expressed as follows: 


r = T [(VJi - (V u ) e ] 
= T (V u ). - T (Vu), 


- Al i - Ar e 


(44) 
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that is, 

4r i - T <T u )i 

*r e - t(t u ). 


The difference in velocity potential on and that on S 2 

at the front stagnation point is equal to Al^ and that at the 
rear stagnation point is equal to Ar e (fig. 22). 

Inasmuch as the flov is considered irrotational, the line 
integral of velocity around any closed path in the flow field must 
equal zero. 



(45) 


where c, for convenience, is taken as any area in the channel 
hounded hy two streamlines and two velocity-potential lines. For 
example. 



The streamlines are so drawn that equal amounts of gas pass 
between them. 



The value of for (i = 1, 2, 3, . . . J; 

where there are j streamlines) is graphically found hy plotting 
pV against - 2 - (fig* 23) and integrating graphically (fig. 24). 
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Then, 


*,=£» (48) 

that is, 

*1 = 0, = £, ^3 = | w for j = 5 

These relations are used to relocate the potential lines and 
the streamlines. The velocities are recomputed and the flow net- 
work again checked. Satisfactory accuracy is usually obtained in 
two adjustments of the flow network. 

The accuracy of the network from 0 4 and the rear stagnation 
point and from 0^ and the front stagnation point cannot be 

checked, inasmuch as the accuracy of the assumed streamline bound- 
ing the channel in these regions cannot be checked by these methods. 
If the network for uniform flow some distance ahead of and behind 
the cascade is, however, smoothly blended with the network between 
the blade surfaces, the surface velocities from 0 4 to the rear 

stagnation point and from 0g to the front stagnation point can 
be approximated with reasonable accuracy. 

The first approximation of the blade shape may be altered by 
increasing the camber or changing the blade thickness or otherwise 
altering the blade shape to satisfy equation (44) and to establish 
a desirable velocity distribution. For example, it is desirable to 
have the maximum velocity on the blade surface no more than 15 per- 
cent greater than the leaving velocity, as specified by the vector 
diagram, inasmuch as too great a reduction in velocity in the 
region of the trailing edge on the suction surface is likely to 
result in flow separation. Accelerating flow over as much of the 
blade surface as possible in order to minimize boundary- layer 
growth is also desirable. 

Blade sections are developed for a number of radii and faired 
to form a complete blade. Stator blades may be designed with a 
linear taper, but the operating stress in a rotor blade for a given 
tip speed will be reduced if a parabolic taper is used. 
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TABLE I - BLADE COORDINATES 



Shaded areas are boundary layer 



Root section 

Pitch 

section 

Tip section 

(in. ) 

x i 

(in.) 

*8 

(in.) 

X 1 

(in. ) 

*2 

(in.) 

x i 

(in. ) 

*2 

(in.) 

0 

0.063 

-0.053 

0.063 

-0*053 

0.063 

-0.053 

.100 

.114 

- .068 

.114 

- .068 

.114 

- .068 

.200 

.148 

- .083 

.148 

- .083 

.148 

- .083 

.300 

.168 

- .110 

.167 

- .105 

.165 

- .100 

.400 

.179 

- .138 

.173 

- .130 

.167 

- .122 

.500 

.172 

- .171 

.164 

- .161 

.156 

- .151 

.600 

.152 

- .211 

.143 

- .1P9 

.134 

- .187 

.700 

.117 

- .258 

.107 

- .234 

.096 

- .229 

.800 

.064 

- .316 

.053 

- .299 

.041 

- .281 

.900 

- .012 

- .388 

- .016 

- .365 

- .029 

- .342 

1.000 

- .116 

- .471 

- .115 

- .442 

- .114 

- .412 

1.100 

- .245 

- .569 

- .229 

- .531 

- .212 

- .493 

1.200 

- .407 

- .686 

- .369 

- .638 

- .331 

- .589 

1.300 

- .612 

- .835 

- .542 

- .769 

- .471 

- .703 

1.400 

- .869 

-1.022 

- .754 

- .930 

- .638 

- .838 

1.500 

1.600 

-1.176 

-1.244 

-1.008 

-1.126 

- .840 
-1.101 

-1.001 

-1.178 


Section 

Pitch 

(in.) 

Chord 
( in. ) 

Axial 
depth 
( in. ) 

Circumferen- 
tial depth 
(in. ) 

Leading- 
edge radius 
(in.) 

Trail ing- 
edge radius 
( in. ) 

Root 

1.193 

2.052 

1.535 

1.275 

0.050 

0.010 

Mean 

1.334 

2.074 

1.585 

1.245 

.050 

.010 

Tip 

1.475 

2.100 

1.638 

1.216 

.050 

.010 
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TABLE II - POSITION OP BLADE PRESSURE TAPS 



Root section 

Pitch 

section 

Tip 

section 

T_ 

X 

X 

v. 

X 

X 

Vn 

X 

X 

n. 

1 

2 

n 

1 

2 

n 

1 

2 

(in. ) 

( in. ) 

( in. ) 

( in. ) 

( in. ) 

( in. ) 

( in. ) 

( in. ) 

( in. ) 

-0.050 

0 

0 

-0.050 

0 

0 

-0.050 

0 

0 

.200 

.150 


.200 

.149 


.200 

.148 


.400 

.178 


.400 

.173 


.400 

.168 


.700 

.119 


.700 

.108 


.700 

.097 


1.000 

-.114 


1.000 

-.115 


1.000 

-.114 


1.300 

-.614 


1.300 

-.543 


•1.300 

-.472 


1.435 

-.970 


1.435 

-.934 


1.435 

-.890 


1.300 


-.835 

1.300 


-.770 

1.300 


-.706 

1.100 


-.570 

1.100 


-.532 

1.100 


-.497 

.800 


-.316 

.800 


-.299 

.800 


-.281 

.500 


-.172 

.500 


-.162 

.500 


-.152 

.300 


-.110 

.300 


-.105 

.300 


-.100 
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* 



Figure 1.- Blade section showing potential lines used 
for calculating the surface velocities. 


Surface critical velocity ratio 
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Blade-surface length, S, in. 

(a) Root section. 


Figure 2.- Blade-surface velocities at design conditions. 
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Surface critical velocity ratio 



2.5 2.0 1.5 1.0 .5 0 .5 1.0 1.5 2.0 2.5 

Blade-surface length, S, in. 

(b) Pitch section. 

Figure 2.- Continued. Blade-surface velocities at design conditions. 
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Figure 2.- Concluded. Blade-3urf ace velocities at design conditions. ® 
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Figure 4.- Survey planes 1*5 chord lengths 
upstream of the turbine blades. (All 
dimensions are given in inches.) 
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Figure 5. - Variation of stagnation and static pressures in the cascade tunnel at 1.5 

chord lengths upstream of the turbine blades. 
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Weight -flow parameter 
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Figure 6. - Variation of weight-flow parameter at cascade entrance 
with ratio of stagnat ion-to-ambient-static pressure across 
cascade • 
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Figure 7. - Measured stagnation and static pressure at 0.1 
chord downstream of cascade. 
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Figure 8. - Comparison of ratios of experimental and design 
stagnation-to-static pressure at 0.1 chord downstream of 
blades. 
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Figure 9. - Variation of gas-discharge angle at 0.1 chord downstream of cascade. 
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Circumferential distance, in* 

(b) Circumferential variation at 
the 10.35-in. radius. 

Figure 9.- Concluded. Variation of 
discharge angle at 0.1 chord down- 
stream of cascade. 
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Figure 10.- Radial variation of critical ratio ^ — and 

v cr v 

tangential component of critical velocity ratio _H_ 
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0.1 chord downstream of cascade. 
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(b) Weight-flow parameter. 

Figure 11.- Radial variation of axial component of velocity 


and weight -flow parameter at 0.1 chord downstream of 
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Figure 12.- Radial equilibrium parameters. 
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Figure 13.- The flow of a gas in the channel 
between curved surfaces. 
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Figure 14.- Chart for finding the average velocity through the blade channel. 
(Two 12-by-20-in. prints of this figure are attached.) 
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Figure 15.- Chart for finding J from channel width and hlade 
curvature. (A 10-hy-l6-in. print of this figure is attached.) 
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Figure l6.- Chart for finding K/J from channel width and "blade 
curvature. (A 10-by-20-in. print of this figure is attached.) 
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C,n 0 /2 

Figure 17.- Variation with channel width and blade curvature of 
ratio of velocities along blade suction surfaces to velocity 
at channel center. (A 10-by-l6-in. print of this figure is 
attached. ) 
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C |ng/2 

Figure 18.- Variation with channel width and blade curvature of 
ratio of velocities along blade pressure surfaces to velocity 
at channel center. (A 10-by-l6-in. print of this figure is 
attached. ) 
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Figure 19.- Assumed camber lines. 
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Figure 20.- Distribution of symmetrical airfoil on assumed camber lines. 
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Figure 21. - Flow network in channel between blades and assumed distribution of weight- 

flow parameter. 
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Figure 22.- Surface velocities and velocity potential. Dashed 
lines are based on an assumed boundary streamline. 
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Figure 23.- Schematic diagram of varia- 
tion of flow intensity with distance 
along velocity-potential line. 
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Figure 24.- Variation in stream function 
with distance along a velocity-potential 
line. 
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(a) Chart for finding JZ m . 
finding average velocity through blade channel. 
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Figure 14. 


Concluded. 
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(b) Chart 
for finding 
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for finding f. 
average velocity 


through blade channel. 
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Figure 15. - Chart for finding J from channel width and blade curvature 
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Figure 16. 
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